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Tests of local hidden variable theories using measurements with continuous variable (CV) out- 
comes are developed, and a comparison of different methods is presented. As examples, we focus 
on multipartite entangled GHZ and cluster states. We suggest a physical process that produces the 
states proposed here, and investigate experiments both with and without binning of the continuous 
variable. In the former case, the Mermin-Klyshko inequalities can be used directly. For unbinned 
outcomes, the moment-based CFRD inequalities are extended to functional inequalities by consid- 
ering arbitrary functions of the measurements at each site. By optimising these functions, we obtain 
more robust violations of local hidden variable theories than with either binning or moments. Recent 
inequalities based on the algebra of quaternions and octonions are compared with these methods. 
Since the prime advantage of CV experiments is to provide a route to highly efficient detection via 
homodyne measurements, we analyse the effect of noise and detection losses in both binned and 
unbinned cases. The CV moment inequalities with an optimal function have greater robustness to 
both loss and noise. This could permit a loophole-free test of Bell inequalities. 



I. INTRODUCTION 

In 1964 Bell showed that the predictions of quantum 
mechanics are not compatible with those of local hid- 
den variable theories (LHV) [l|. Since then, there have 
been experimental tests [2|-lj| of these predictions, which 
can be used to rule out some LHV alternatives to quan- 
tum mechanics (QM). While experiments to date support 
quantum mechanics, the experimental strategies used 
have loopholes [5|. In some experiments, problems are 
caused by insufficient spatial separation of measurements 
— so that the assumption of local causality that under- 
pins the LHV model cannot be applied to start with. In 
other experiments, the poor photon detection efficiencies 
[7| mean that correlations are not strong enough to vio- 
late a Bell inequality, unless additional assumptions are 
utilized. In short, a clear violation applicable to all local 
hidden variable theories has not yet been experimentally 
observed. 

It has been widely suggested that the use of highly ef- 
ficient detection techniques like homodyne detection |8| 
may provide an avenue to overcome the detection effi- 
ciency problems [^ with Bell's theorem. Already, contin- 
uous variable (CV) realizations of the Einstein-Podolsky- 
Rosen (EPR) paradox have been achieved with high ef- 
ficiencies y, using optical homodyne detection [10|. In 
this case, the EPR paradox is reformulated in terms of 
field quadrature phase amplitudes which have commuta- 
tors like position and momentum. The EPR paradox, 
however, only demonstrates the incompatibility of local 
realism (LR) with the completeness of quantum mechan- 
ics. It does not test local realism itself. Proposals to 
demonstrate an incompatibility of quantum theory with 
local hidden variable theories using continuous variables 
require new formalisms, since Bell's original inequality 
was only applicable to discrete or binned outcomes. 

Early proposals for such tests using CV measurements 
[8| relied on binning the initially continuous outcome do- 



mains to obtain binary results. Here the original Bell 
Clauser-Horne-Shimony-Holt inequalities |ll| are appli- 
cable. However, these predicted violations of Bell's in- 
equality were small and involved states that were diffi- 
cult to prepare. Wenger et al. advocated a root bin- 
ning technique to enhance the level of violation |13 |. 
Nha and Carmichael [lj| proposed a conditional state 
preparation to obtain the required nonclassical states. 
Munro et al. [1J| and Acin et al. [l^ suggested us- 
ing binned CV measurements made on N spatially sep- 
arated qubit systems prep ared in a Greenberger-Horne- 
Zeilinger (GHZ) state [1^. and noted improved violation 
of the Mermin-Ardehali-Belinskii-Klyshko (MABK) Bell 
inequalities [ItI as N is increased. 

More recently^Cavalcanti, Foster, Reid and Drum- 
mond (CFRD) [l8| have derived a continuous variable 
inequality for LHV theories, formulated directly in terms 
of moments, without binning. This allows a simple test of 
local hidden variable theories for genuine CV outcomes. 
Violations of this CV moment inequality were found pos- 
sible for GHZ states involving a large number of sites 
(A^ > 10). Based on the CFRD inequalities, Shchukin 
and Vogel |19| introduced further inequalities considering 
4 or 8 observables at each site. Salles et al have shown^] 
that at least three sites are necessary for the CFRD in- 
equalities to be violated, and have more recently gen- 
eralized this technique to multilinear mappings [21|. In 
general, the use of ordinary moments has the advantage 
that it is more robust against loss and noise than earlier 
binning or conditional techniques. However, the large 
number of multipartite sites involved presents practical 
difficulties. 

In this paper, we develop functional CFRD Bell in- 
equalities for testing local hidden variable theories, us- 
ing arbitrary functions of continuous variable outcomes 
at each site. Variational calculus is used to optimize the 
choice of measured function, to maximise the violation of 
the Bell inequality. As a result, we obtain predictions for 



violations of LHV models with fewer sites, larger losses 
and greater degrees of decoherence than previously pos- 
sible. This approach of optimising the function of a mea- 
surement result for a particular purpose should be gener- 
ally useful in the field of quantum information. Prelim- 
inary results have been presented as a short Letter |23 |. 
Here we provide a more detailed and thorough account 
of this procedure. 

The work of Acin et al. [I^ and Cavalcanti et al. [l8fl 
has revealed an exponential improvement in Bell inequal- 
ity violations as the number of sites N increases, for CV 
measurements. These results support the earlier work 
of MABK, which also showed an exponential improve- 
ment of violation for iV-partite states for spin measure- 
ments. Since these inequalities require detection at every 
site, one might expect this improvement to be counter- 
acted by increasing sensitivity to detection losses. Work 
on the robustness of Mermin-Bell and CFRD-Bell viola- 
tions with loss [ISl . |24| shows that this is not the case. 
Other calculations using different LHV assumptions (25| 
have suggested a loss threshold as low as "qcrit = 0.5 as 
N ^ oo. However, since this proposal uses additional 
assumptions about the LHV models, it may not rule out 
all LHV theories. 

In the present paper, we introduce a new approach to 
nonlocality that utilizes a general functional optimiza- 
tion of continuous variable observables. We find the 
optimal function that maximises a violation of the CV 
Bell inequality for a given efhciency 77 and state purity 
p. We show that the optimal function has the form 
x/{l + Snx^), where e^r is a parameter related to N and 
77. This gives an inequality for all LHV models, which is 
violated by the GHZ states of (H^]) for A^ > 5. The viola- 
tion increases exponentially with N . For a pure state the 
loss threshold rjcrit decreases asymptotically to 7700 =0.69 
as N ^ 00, in the most symmetric case. Thus both the 
number of modes required, and the required efficiency are 
dramatically reduced. The wider range of quantum states 
considered here include states which can be straightfor- 
wardly obtained from a polarisation GHZ state of the 
type generated in current experiments. These are dif- 
ferent to the extreme photon-number correlated states 
considered by Acin et al, which may be an experimental 
advantage. 

A simple example of a functional transformation is ob- 
tained when binning of a CV observable is used to give 
discrete binary outcomes Q. In this case, the CFRD 
inequalities reduce to those of Mermin J16ll . Here, our 
analysis is similar to that of Acin et al. |15| . For binary 
outcomes, one can use extensions of Mermin's inequali- 
ties [131 for the extreme photon-number correlated states 
|7V) = (|0)®^ + |1)®^) /V2. One then obtains an expo- 
nential increase of violation of LHV with N, and violation 
for A^ > 3. In addition, we extend the analysis of Mermin 
and Acin et al. [15|, and calculate results for homodyne 
detection for more feasible types of state. However, we 



using functional moment inequalities. The critical effi- 
ciency for a pure state at large N is r] = 0.79, which 
is substantially higher than for the optimal functional 
CFRD approach. 

The organisation of the paper is as follows. Section 
II gives the general notation used for continuous vari- 
able moments, and reviews the original CFRD inequality. 
In Section HI we review and extend the CFRD inequal- 
ity for the case of quadrature measurements on general- 
ized GHZ states. Section IV considers the multipartite 
binning approach that uses the MABK Bell inequalities, 
which was also studied recently by Acin et al. for quadra- 
ture detection. In both Sections HI and IV, we consider 
several types of GHZ state, and also present results for 
cluster and W-states [26l. [27|. Next, in Section V these 
approaches are generalized to include functions of the 
observables. We compute the maximum violation of the 
inequality for several choices of function, in the case of a 
generalized multipartite GHZ state. Functional transfor- 
mations using an optimum function, i.e., x/{l + Enx^), 
are shown to give a substantial improvement to results 
obtained with the CFRD approach. In Section VI these 
results are compared with the Shchukin-Vogel (SV) in- 
equality for larger numbers of observables, showing that 
this gives identical results to the simpler CFRD method, 
in the cases studied. In Section VII we include the ef- 
fects of loss and noise on the violations of the functional 
moment inequalities and MABK inequalities. Optimiz- 
ing the choice of measured function reveals a quantum 
nonlocality for larger losses and for greater degrees of de- 
coherence than possible previously. Section VIII gives an 
overview of the results. 

In summary, we examine in this paper the effect of 
loss and state preparation impurity on tests of local re- 
alism using functional moments. Our motivation is the 
fundamental challenge of developing an experimentally 
feasible loop-hole free test of local realism, which bal- 
ances the need for a simple experimental realization with 
robustness against loss and noise. 
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proach is more susceptible to loss and decoherence than 



II. CONTINUOUS VARIABLE 
MEASUREMENTS 



Early work on Bell inequalities utilized measurements 
of spin, making it seem that obtaining discrete measure- 
ment outcomes was a crucial ingredient to violating a 
Bell inequality. In fact, this is not essential. The CFRD 
inequality is a multipartite test of local realism involving 
ANY two observables measured at each of N sites, with 
causal separation between measurement at each site. 



A. General CFRD inequality 

It was shown in [18| that the following inequality must 
hold for all LHV theories: 
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(2.1) 



where Xi,., Pn are the outcomes of two arbitrary measure- 
ments |23| at each site labeled n. The measurements are 
causally separated, and ()p indicates an average over a 
set of local hidden variables A with probability P (A). 

In quantum mechanics, these correlations involve op- 
erators a;„, Pn at each site n, with eigenvalues x„, p„. 
Unlike in hidden variable theories, these inequalities can 
be violated for certain quantum states, if the operators 
do not commute. Despite the notation used, nothing is 
assumed about the type of measurement operators. In 
LHV theories, the inequality must be satisfied by mea- 
surements with any type of spectrum, discrete or contin- 
uous, and even unbounded outcomes. For this inequality 
to be useful, one must find states and observables that 
violate the inequality within quantum theory, and are 
obtainable in an experiment. 

It is important to recall here that when x„, p„ are 
described within quantum mechanics by non-commuting 
operators, the corresponding physical variables are not 
simultaneously measurable. Thus, the TV-fold products in 
the above inequality each involve a sum over 2^ distinct 
correlation measurements with different settings of the 
measurement apparata. 



units, (see the Appendix) the position and momentum of 
a harmonic oscillator, with comnmtators: 
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One can adjust the local oscillator phase to detect linear 
combinations of position and momentum. Hence, we can 
introduce a general quadrature phase operator: 



X„ = Xn cos 9 + Pn siu 6 
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For continuous measurements of quadrature variables 
on entangled boson states, we define generalized GHZ 
states to be: 
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Here the integer r is the maximum number of modes 
having nonzero photon number. The critical efficiency 
Vcrit required for violation of the CFRD inequality in the 
symmetric case of r = N/2 tends to r/crit = 0.81, as 
A^ ^- oo. Quadrature measurements with local oscilla- 
tors are inherently highly efficient compared to spin or 
photon-number measurements, with reported efficiencies 
of 99%. However, generation losses from mode-matching 
can degrade the experimental efficiency, so 81% is still 
a challenging practical benchmark, especially with large 
numbers of correlated modes. 



Functional transformations 



B. Quadrature measurements 

While the inequality is valid in general, we wish to ap- 
ply it to experiments in quantum or atom optics, which 
are known to be able to produce highly correlated, non- 
classical quantum states. Operationally, we will consider 
the case where conventional photodetectors are replaced 
by a homodyne detector using a local oscillator, in order 
to detect field quadratures rather than photon numbers. 
As well as involving an intrinsically efficient experimen- 
tal technique, this method is already known to be robust 
against loss, without additional assumptions. Although 
quadrature measurements are traditionally used for co- 
herent and squeezed states, they of course can also be 
used for states with well-defined overall particle number. 
In principle, exactly the same techniques apply to any 
boson fields, including, for example, ultra-cold atoms. 

We assume there are A^ sites at which one can make 
observations, with an adjustable local oscillator phase 
to define local quadrature operators x„, p„. In the 
case that Xn and p„ are associated with single-mode 
bosonic creation and annihilation operators aj^, a„, one 
has Xn ~ (a„ + a,\) /2 and p„ = i{al^ — a„)/2. These 
are the quadrature phase amplitudes associated with the 
boson annihilation operator a„, or for suitable choice of 



Functional transformations of any output measure- 
ments are always possible experimentally. Any data file 
of experimental results is trivially processed into a func- 
tion or set of functions of the original measured data. 
While this procedure is not very useful for binary results, 
it represents a large functional space of measurements in 
the continuous variable case. 

In terms of hidden- variable theories, this simply means 
that all measured quantities x can be replaced by a func- 
tion f{x). In quantum mechanics, this implies that all 
eigenvalues e are replaced by /(e), or in operator lan- 
guage, the measurement operator x is replaced by / {x) . 
The details of this procedure are explained in later sec- 
tions. We will show that such functional transformations 
greatly reduce the number of sites and the critical effi- 
ciencies needed to violate a Bell inequality. 



III. CFRD INEQUALITY 

In this section we will review the CFRD inequality |18fl 
as applied to quadrature measurements, and provide a 
more complete analysis of its violation than in the Letter 
that introduced it. Quadrature observables at A^ sites 
are considered here, as described in the previous section. 



A. Quadrature CFRD inequality 

At each site, the sign associated with p„ can be chosen 
to obtain either a;„ + ip^ = On or x„ — ip„ = aj^ for each 
of the terms in the left hand side {LHS) of (P7T|) . The 
choice of sign is important to maximise the violations for 
a given state. Defining a variable s„ S {—1,1} to repre- 
sent this choice for each site, we denote by A„(— 1) = a„ 
and An{l) = aj^ the corresponding combination of mea- 
surements at each site, 
duces to 



The CFRD inequality then re- 
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CFRD showed that within quantum theory, the symmet- 
ric state 

violates this inequality for N > 10. Note here that there 
is a spatial mode index that has been dropped for con- 
venience but is given by the order of the kets, so that 
this state represents A^ distinct modes at spatially sep- 
arated sites, and each mode is occupied by one or zero 
photons. Thus, the states |0), |1) are eigenstates of the 
number operator a^a with corresponding eigenvalues 
and 1 respectively, so the prediction could in principle 
be tested with photonic GHZ states produced in the lab- 
oratory. While a value of iV = 10 is not impossible, it 
is a large number of correlated modes for a practical ex- 
periment which involves multiple down-conversions and 
beam-splitters. 

This type of symmetric GHZ state is obtainable from 
the usual GHZ polarized state via a local unitary trans- 
formation at each location. For example, the symmetric 
state given above can be prepared from an iV/2-photon 
polarization GHZ state, of a type that is found in some 
current experiments |28| - |3C| 



\GHZP) = -^ (|H)^^/2 + \Vf^'"' 
V 2 ^ 



(3.3) 



A possible procedure is shown in Fig. [TJ where |7J), \V) 
represent horizontally or vertically polarized single- 
photon states, by passing each photon through a po- 
larizing beam splitter (PBS). This converts polarization 
modes into spatial modes, so that \H) = \1)h\Q)v and 
\V) = |0)_f/|l)v where \Q)h/v and \1)h/v are number 
states for the horizontally and vertically polarised modes 
respectively. The \GHZP) state then becomes |-y) of 

(I32D. 

On-demand state generation 

Some caution is needed here in interpreting current 
GHZ experiments [23 ■ These typically involve condi- 
tional measurements, rather than a known initial quan- 
tum state. To our knowledge, current experiments only 
generate GHZ states in the sense that state is identified 



after some number of photon-counting events have al- 
ready taken place. Here we assume a more traditional 
quantum mechanical state preparation, in that we as- 
sume a GHZ state can be generated on demand. We 
note that a proposal already exists for carrying out on- 
demand generation of a singlet state [31j, which is the 
simplest case. Current experimental techniques for GHZ 
states would therefore need similar modifications to pro- 
duce a well-defined GHZ state 'on demand'. 
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Figure 1: (Color online) Schematic of an experimental setup 
to test the CFRD inequalities. We assume an A'^/2-photon 
GHZ state can be generated. Then by using a PBS, each 
photon is sent to one of the two output spatial modes de- 
pending on its polarisation. The entanglement in polarization 
is therefore translated into entanglement in spatial modes to 
produce the symmetric states \N/2) with A'^ field modes. At 
each site, high-efficiency homodyne detection is used to mea- 
sure the correlations. 



B. Generalized GHZ states 

A class of states that violate the CFRD inequalities is 
represented by the entangled states 

IV;) = ci|0)«^ni)^^^"''^ + C2|l)^''|0)^(^-'') . (3.4) 

Thus r = N corresponds to extreme photon-number cor- 
related states, a superposition of a state with photons 
at all sites and a state with 1 photon at each site, while 
other values of r correspond to intermediate cases. These 
can of course be transformed into each other with lo- 
cal unitary operations, although in practical terms such 
transformations are not always feasible without losses. 

Within quantum mechanics, the LHS of the CFRD 
inequality ()3.ip is non-zero with the choice 



LHS = ( a] ..alar-\-i...aN 
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Here, |cip + |c2p = 1, and the violation is maximized for 
ci = C2 = l/\/2. For the right hand side (RHS) of the 



inequality, we obtain 

RHS = ((a|ai + l/2)...(aJYaAr + l/2)) 

= c\{l/2Y{i/2f-'- + c2(l/2)^-''(3/2)'^(3.6) 
If we use photon-number correlated states with r = N/2, 



RHS = {\r'\ 



(3.7) 



which is independent of the amplitudes ci , C2 . 

Alternately, if we fix ci — C2 — l/v2 but change r, 
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There is a violation of the inequality when the Bell 
observable Bn = LHS/RHS > 1. We note from Fig. 
[5] that no violations are possible for r = N, and the 
optimal case has r approximately equal to N/2. For N 
even, where r = N/2 there is a violation if 
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in agreement with Cavalcanti et al. [18| who have shown 
this requires N > 10. For odd N, we select r = {N±l)/2 
to obtain a violation if 



2JV-I ^ '^(N^l)/2 _(_ 3(JV+l)/2 

which is obtained for iV > 11. 



(3.10) 




Figure 2: (Color online) Bell violation as a function of r with 
ci = C2 — I/V2. Different A^ (10-14) values are distinguished 
by their line style. The allowed values of r are of course only 
the natural numbers; the lines are drawn to guide the eye. 



either because they are simpler to generate experimen- 
tally, or are more resistant to certain forms of decoher- 
ence. In addition to the GHZ states [2^, four and six- 
photon cluster states have been realized experimentally 
|29l ISOJ . We thus present initial results for the CFRD 
inequality using cluster and W-states. 

Consider first the four-qubit cluster state [23, [2^, |30| 
defined as 



10 = -(|1)|1)|1)|1) + |1)|1)|0)|0) 

+ |1)|1)|0)|0)-|0)|0)|0)|0)) (3.11) 

and also consider the four-qubit W-state [2g 



m = -(|1)|0)|0)|0)-I-|0)|1)|0)|0) 

+ |0)|0)|1)|0) + |0)|0)|0)|1)). (3.12) 

In this case the qubit values correspond to the occupa- 
tion number (0 or 1) of a spatial mode, as described in 
Sections III A and IIIB. The cluster state of (|3.1ip is a su- 
perposition of two terms which correspond to GHZ states 
of type p.4p with r = N/2 and r — N respectively. From 
Fig. [21 we note there is no Bell violation when r = N. 
The effect is that this cluster state case is less optimal 
than the symmetric GHZ state (r = N/2) and no vio- 
lation of the iV = 4 CFRD inequality is observed. For 
W-states ((3T2| . the value of the LHS of (pTS]) is always 
so that no violations of this particular CFRD inequality 
are found. Similarly negative results were obtained for 
the arbitrary N extensions of these superpositions lip) — 

where N was increased. 

We can also consider the specific four-photon clus- 
ter state that has been produced experimentally by P. 
Walther et al El: 



\C) = -i\H)\H)\H)\H) + \H)\H)\V)\V) 

+ \V)\V)\H)\H)-\V)\V)\V)\V)), (3.13) 

where \H) and \V) represent respectively horizontally 
and vertically polarized single-photon states. Thus, 
\H) = \\)h\^)v and \V) = \Q)h\1)v, the \Q)h,v and 
\1)h/v being number states for each polarized mode de- 
noted by subscript H and V . The state (I3.17P is a four 
photon state and the qubit value is encoded in the polar- 
isation. We also consider the four-photon W-state: 



C. Cluster and W states 

Other forms of multipartite entangled states, such as 
W-states and cluster states [2a, u3, are also of interest. 



\W) = -{\V)\H)\H)\H) 
+ \H)\H)\V)\H) 



^\H)\V)\H)\H) 
\H)\H)\H)\V)).{3.U) 



As explained in Section 3A, these states can be converted 
to the states |^) of (|3.2p with N = 8 via the apparatus 



of Fig. 1. We can then use the N=8 CFRD inequahty 
for testing the nonlocal behaviour of these states, which 
are written explicitly as 
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\W) = -(|0)|1)|1)|0)|1)|0)|1)|0) + |1)|0)|0)|1)|1)|0)|1)|0) 



-|1)|0)|1)|0)|0)|1)|1)|0) 
-|1>|0)|1)|0)|1)|0)|0)|1)). 



(3.16) 



Evaluation of the left- and right- hand sides of the CFRD 
Bell inequahty (gH) gives: LHS = 1/16 and RHS = 
3^/2^. Thus, no violations are possible for the N = 8 
CFRD inequality with these 4 photon states so that in 
this case the cluster state is less optimal than the r — N/2 
GHZ state (|3.4p . However if we consider the N photon 
extension of the 4-photon states (|3.17l) and (13. 15^ : 



\C) 



1 






we find a more positive result: LHS = 1/16 and RHS = 
(3/4)^/^. so that a violation is possible when iV > 20. 
We note in this case no violations are obtained for the 
W states. 



D. Angular dependence 

The CFRD inequality holds for arbitrary measure- 
ments and is therefore not restricted to the orthogonal 
choice specified in the previous section. Noting that Bell 
inequalities for correlated spins are optimized at polar- 
izer angles that are distinct at different sites, we wish 
to investigate the effect of an arbitrary local oscillator 
phase. 

Using general quadrature operators X^^ = (a„e^*^" + 

aJje*^")/2 for each site n, and defining F„ — Xf^" +iXn^, 
we substitute these in the CFRD inequality (|2.ip with 
Xn = X^ and p„ = X^ . The inequality then becomes 
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which gives additional terms in the RHS. We note that 
in QM, these extra terms correspond to operators that 



change the boson number by two. Hence for states like 
|0011) -t- 1 1100) they will be zero, leaving the RHS in- 
variant with the choice of phases On , 0'^ . Therefore, for 
GHZ states, the RHS will be just as in Eq ((XB)) . 

For simplicity, we treat the case of iV = 2 as an il- 
lustration. Within QM, the LHS for the Beh state 
|V)=ci|G)|l)+C2|l)|0)is: 

LHS = \{{X<l + tXf){xt + ^xf)\' 

= \{[{xi cos6 + pi sin 6') + i{xi cos9' -f pi sin6'')] 
X [(2:2 cost/) + P2sin0) + i{x2COS(j)' +p2 sin (/>')]) |^ 



\-ciC2[icos{6 - 4>) - cos 
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since here (xia;2) = {P1P2) = C1C2/2 and {X1P2) = 
{p2Xi) = 0. The LHS is maximised either with the 
orthogonal choice 9 = (f) = 0, 6' = — 7r/2 = —(f>', for 

which (Im{A-F2}) = 0, (Re{Fii^2}) = C1C2; or with the 
choice 9 = 0, 9' ^ 7r/2, <j) = 7r/4, <j)' = -7r/4 for which 

(Im{Ai^2}) = C1C2/V2, (Re{AA}) = C1C2/V2. In both 
cases, the LHS is equal to |ciC2p, as in Eq (|3.5p . Hence, 
just as in the orthogonal quadrature case, there is no 
violation. 

The case for the extreme photon-number correlated 
state \ijj) = Ci|0)|0) + C2|l)|l) is even less favourable. 
The value of LHS can still be maximized to |ciC2p, with 
angle choice 9 = (/) = 0, 9' = cj)' = ±tt/2 ot 9 = 0, 
9' = 7r/2, — — 7r/4, (f)' — tt/4. However, the value of 
the RHS is increased even further, so that again, there 
is no violation. 

More generally, the inequality is violated for N > 10, 
that is, with more correlated modes and observers. We 
find that the optimal violations are found with phase 
choices: 6*1 = • • • = 0„ = 0, and 9^ — ■ ■ ■ — 9^ — —9^,^^ — 

• • • = -9'n = -^/2 or 9n = (-l)"+i7r(n - 1)/(27V), 9^ = 
9n + tt/2 {n < r) and 9n = {-l)"n{n - 1)/(27V), 9',, = 
9n - tt/2 {n>r). The LHS and RHS then become, for 
arbitrary TV and r, with ci = C2 = \/ \f2: 



RH S max 



I C1C2 
oN-r 



I'- 1/4 , 
+ 3'^ 



2^+1 



(3.20) 



In summary, we find that the violations are indepen- 
dent of the relative phases between different sites [22| — 
unlike the usual case with Bell inequalities for correlated 
spins — , thus confirming a result pointed out in [l8| for 
this larger range of states. 



IV. MABK INEQUALITIES 

In early proposals for obtaining inequalities to test lo- 
cal realism with continuous variable measurements, only 
the sign of the quadrature variable is recorded Q . As far 



as LHV theories are concerned, it is just the final out- 
comes that are relevant, not the quantum operators cor- 
responding to them. This allows one to use the same logic 
that is employed for dichotomic spin variables. Our re- 
sults here are therefore similar to those of earlier workers 
B H, llH ■ As in the work of [1^ , we find an exponential 
improvement in Bell violation with increasing number of 
sites, but we demonstrate that result for a larger class of 
states, some of which may be within experimental feasi- 
bility. 



, for A^ even, and 

Sn - 2-(^-i)/2Re(Im) {{nN)p} (4.4) 

for N odd. The inequality also holds for Sn — 
2-(^-i)/2 ^^(Re{nAr})2 + (Im{njv})2^ . These in- 
equalities have been considered recently for the case of 
GHZ states ^^, where r = TV, by Acin et al [Hi. We 
can also define Hn with the adjoint operators at any of 
the sites, to obtain the same inequality. 



MABK approach with binned variables 



B. Three-site example 



Violations of the MABK inequalities, which were orig- 
inally used for spin variables, can also be obtained by us- 
ing binned outcomes of continuous variable observables. 
These are applicable to any quantum state. Here the 
binned measurement x^'" = 20 (.t„) — 1 is the observ- 



able with outcome 



= +1, if Xn > 0, and x„ 



bin 



1, 



if Xn < 0. We define p^"'- similarly, together with cor- 
responding quantum operators x^^ — 20 (i;„) — 1 and 
pbi7i _ 20 (j5^) _ I Next, we can introduce a new 
binned complex observable a^'" at location n through its 
real and imaginary parts, so that Re{a^j*"} = x^™ and 
ImjaJ^™} = pj^"'. There is also a corresponding operator. 



xhin 



^bin 



For N sites, we can modify the CFRD proof directly, 
assuming local hidden variables, to get: 



TV 



^bin bin 



•■)pr< n(K"T + [P"l ) 



-.N 



\n=l 



This follows because ([a 



6ml ^ 



/ \ r^bir. 



= 1. 



(4.1) 
The in- 



equality remains valid if any of the complex variables a^"* 
are replaced by their conjugates. 

The CFRD inequality for this binned case then reduces 
to that of Mermin and Ardehali [1^1 by taking the square 
root of both of sides. However, the inequality was later 
strengthened by Belinskii and KlyshkofT^, to give an ap- 
proach which we refer to as the MABK inequality. This 
approach can be extended to account for more general 
angles and states, by introducing binned angular observ- 
ables Xf,'''™ = 20 iXf) - 1. Next, a^*" is replaced with 



F°™ = X^'''™ + iXl '''™. The notation 9^ and 0'^ repre- 
sents the angles defining the two local oscillator phases 
chosen at site n, where k = l,...,iV, giving 2^ distinct 



sets of measurement settings. 
Then, defining Hat = f^'"F^'". 



pbin 



we can test for 



violations of the even stronger MABK inequality 



Within quantum mechanics, the LHS of the Bell in- 
equality (|4.2p can be calculated via binned quadrature 
phase amplitudes. It is instructive to give an example of 
the calculation in detail. For the state (12. 3p . the lowest 
values of A^ and r that give a violation are A^ = 3 and 
r = 1: 

|V') = ^(|0)i|l)2|l)3 + |l>l|0)2|0)3). (4.5) 

Using the harmonic oscillator wavefunctions given in the 
Appendix, we find that: 

(Hg) = (^pbinpMnpMn'^ 

= (-)3/2[cos(e' - - 7) - COs(e -(f)' ^ i) 

-cos{e' - 0' - 7) - cos{e' - - 7') 
+i{cos{e - 0' - 7) + cos[e - - 7') 
+cos{e' - (/) - 7) - cos{e' - 0' - 7'))] ,(4.6) 

where we denote 9i — 6. 62 — (j), and ^3 = 7. 

Choosing the two measurements made at each site to 
be orthogonal, and measurements at all sites along the 
same (or opposite) directions, so that 6* = (/> = 7 = 0, 

and 9' = —((>' = —7' = — 7r/2 (which means that H3 = 

alinhmCbtn), wc obtain 



^3 = 2-(3-l)/2^(Rc{H3})2+0 = 2(2/^)3/2 > 1 



(4.7) 



Alternatively, using different measurement angles: 9 — 

0, = -7r/6, 7 = -27r/6, 9' = 7r/2, 0' = -47r/6, 7' = 
— 57r/6, we can rotate the LHS moment, so that 

^3 = 2-(3-i)/2^o+(Im{H3})2 = 2(2/^)3/2 > 1 . 

(4.8) 

Next, consider the extreme GHZ state of form (|2.3p 
with r = 3: 



where 



'Af 



-N/2 



5jv| < 1 , (4.2) 



(Re{HA,}±Im{HAr})p (4.3) 



|V') = ^(|0)i|0)2|0)3 + |l)l|l)2|l)3). (4.9) 

We obtain a new (H3) with cos (0-1- 0-1- 7), cos{9 + (f)' + j'), 
etc. replacing cos{9 — (j) — 7), cos{9 — (j)' — 7'), etc. in 



8 



Eq. (|4.6p . In this case the choice 9 = 0, 4> = 7r/6, 7 = 
27r/6, 0' = 7r/2, </()' = 47r/6, 7' = 57r/6, gives the maxi- 
mum violation of the inequahty, 



^3 = 2-(3-l)/V(Rc{n3})2 + = 2(2/^)3/2 > 1 , 

(4.10) 
as has been presented by Acin et al. 



C. N-site inequalities 

Continuing in this fashion for higher TV, we find that 

sn = ^{-r'' (4.11) 



for the choice 6'„ = (-l)^+i7r(fc - 1)/(27V), 61^, == 6i„ + 
tt/2 for A: < r, and 9,, = (-l)^7r(n - l)/(2iV), 9'^ = 
9n — ''■/2 for n > r. There is a violation of the Bell 
inequality when Sn > 1. This result has been presented 
by Acin et al. for the case of r = N. We confirm the 
exponential increase with number of sites N, but also 
make the observation that the violation occurs for all 
types of states of the form (|2.3|) . independently of r. This 
contrasts with the result for the CFRD inequality, which 
is strongly dependent on r, requiring r ~ A^/2. Violation 
of the BI with binning is therefore possible in principle 
for A^ > 3, but, as shown in later sections, this strategy 
is very sensitive to losses and noise. 



be examined. We leave open the question of whether other 
Bell inequalities might be more sensitive for the cluster 
and W-states with continuous variables measurements, 
as shown by Scarani et al [32| for spin measurements. 



V. FUNCTIONAL OPTIMIZATION 

In this section, we show that greatly improved experi- 
mental simplicity, together with robustness against deco- 
herence is possible by testing local realism using a func- 
tional moment inequality approach. The advantage of 
the CFRD moment approach compared to the binning 
approach [8| is that it gives a low weight to quadrature 
signals with low amplitude, which reduces the sensitivity 
to noise. On the other hand, binning has the advan- 
tage that it saturates, which increases the relative size 
of the violations. We show that an approach that com- 
bines these two features using functions of quadrature 
variables is possible. 

We analytically calculate the optimal function using a 
variational calculus method. This produces an inequality 
which is violated by the states of Eq. p.4p for A^ > 5, 
i.e. half the number of modes required with the method 
of [1^. The violation increases exponentially with N, 
and we will see that the critical detection efficiency r]crit 
decreases asymptotically to 0.69, a significant reduction. 



A. The functional inequalities 



D. Cluster and W states 

We also present initial results for the violations of the 
MABK CV inequalities using the specific cluster and W- 
states (13. lip and (|3.12p We first consider the case N = 
4 with CV binned outcomes. Using the same method 
as above, and selecting the optimal choice of angles, we 
obtain the best N=4 MABK correlations of 



|^4-C| 
IS'4-wl 



^ V2/4 

4 Vtt 
= < \S4-GHZ 



<\Si- 



i-GHZ\ 



2 \7T 



for cluster and W-states respectively, and thus we see no 
violation here for A = 4. However, on examining the 
generalisation 



Functional inequalities are already implicit in the orig- 
inal derivation of CFRD [1^. For completeness, we 
present a proof that takes explicit account of functions 
of measurements that can be made at each of A^ causally 
separated sites. We consider the measurable joint prob- 



ability F(Af,A|...) for outcome Af, A| 



at loca- 



tions 1,2,.. ., respectively, where 9, (p, ... represents a 
choice of measurement parameter. For local hidden vari- 
able (LHV) theories, the joint probability is written [l| 
in terms of variables A, as 



1 

(4.12) P{xlxt,..:) = j^ 



d\P{X)P{Xl\\)P{xt\\)... , (5.1) 



IV') = (|i 



(4.13) 
of the superpositions p. lip and (|3.12p for higher A^, we 

find |S'jv-c| — 4 (7) so that violations are possible 
for A^ > 10. Calculation for the W-states shows no vio- 
lation for any A^. These calculations are specific to the 
case where the qubit is realized as a photon number of 
1 or at a given spatial mode. The cases p.l7p - p.l6p 
where the qubit is realized in terms of horizontal or ver- 
tical polarisation, as explained in Section III, could also 



where P{Xf\\) is the probability for result Xf given that 
the system is specified by the hidden variable state A. We 
note that while we use the term "hidden variable state", 
the states can in fact be quantum states, denoted p\, 
as long as they are separable: i.e., not entangled. The 
assumption of local causality (or separability) allows the 
factorization in the integrand. From this, it also follows 
that we can write the measurable moments (A"iA2A3...), 
where A„ is an observable associated with a measurement 
at site 71, etc., as 



(AfX|X 



4>^-t 



jP= I d\P{\){X<>)^{Xt)x{X2)x. 



(5.2) 



Here {X^x is the average of X^ given the local hidden 
variable state A, which means that 



{xi) 



dxip{xi\x)x: 



(5.3) 



where X^ denotes the outcome of measurement at site n 
with phase 0. 

Next we construct, for each site n, real functions of the 
two observables fn{X^), gn{X^^ ), and define the complex 
function: 



F„ 



UX'J+ig„{X^;). 



(5.4) 



The complex moment {FiF2...F]\[) can be expressed 
using sums of real-valued expressions of the type 

{h{Xf)g2{X$')...fN{X'^)). Of course, /„(X^) is just 
another observable composed of the observable X^ plus 
some local post-measurement processing. This post- 
processing can be carried out at any time subsequent to 
the measurement. Eq. (|5.2p is therefore valid for 



ifiXf) . . . f{X%))p = J^ rfAP(A)(/(X«))A . . . {f{X^))x. 

(5.5) 
The expectation value of products of the Fn must satisfy: 



{F^...Fn)p= j dXP{X){F,)x...{FN)x, (5.6) 

where {F„)x = (/„(^fJ>A + i{gn{X^))x- Hence, from 
(|5.6p . the following inequality holds: 



\{FiF2...Fn)p\^ < f dXP{X)\{Fi)x\''...\{FN)xfi5.7) 

For any particular value of A, the statistics predicted 
for fn{Xn) must have a non-negative variance, i.e., 
{fn{Xn))l < {.fniXnf)x. Then wc cau rewrite dEZ]) 
explicitly in terms of the /„'s. Using this variance in- 
equality we arrive at a functional moment inequality: 



B 



rN 



n::=i /„(x,^)+^.9„(x«) 



tN 



uLAux^r+9n{x?:) 



< 1 



(5.8) 



We will measure the violation of this inequality by the 
Bell observable B defined above, so that failure of LHV 
is demonstrated when B > 1. 



B. Functional inequality examples 

First, we investigate three specific types of function: 

1. Fractional order moments: fn{x) — gn{x) - 
\x\^'"sign{x). 



I 



DQ 



Triple Binning (s=0.11) 
■ Tanli function (q=4) 
Binned outcomes (s=0,or q»1) 




Figure 3: (Color online) Maximum violations of MABK Bell 
inequalities as a function of the number of modes A'' with 
ideal detection efficiency. Violation of the Bell inequalities is 
achieved when B > 1. The triple step function case is the 
dotted line with s = 0.11, which reduces to the binned case 
when s = (solid line); the tanh function case is the dashed 
line with q — 4, m = 1.4, which reduces to the binned case 
when q is large. Here the parameters m, q, and s are the 
optimal choices to maximize the Bell value. 



2. Triple binning: 
function equal to 
and elsewhere. 



fn{x) is the triple- valued step 
— 1 when X < — s, +1 when x > s, 



3. Powers of tanh functions: 
sign{x)\ tanh{qx)\™ . 



fn{x) = gn{x) 



For all of these functions, we obtain an LHV violation 
when iV > 5, with optimal choice of rn, s, and q. This is 
a much more feasible experiment than that for an integer 
moment with m = 1, which required iV > 10 in CFRD. 
The reason for this is that these three functions all com- 
bine high intensity saturation with small-noise insensi- 
tivity. The region close to is relatively flat, rather than 
discontinuous as in the binning approach. This leads to 
reduced sensitivity to small noise effects, which otherwise 
will lead to a greatly amplified contribution of random 
vacuum noise to the correlations. Here, we again select 
the optimal r at r = N/2 (TV is even) or at r = {N±l)/2 
{N is odd). The extreme state where r = iV is less opti- 
mal. 

It is worth noting that bounded CV observables (cases 
2 and 3 above) can also be used with the MABK in- 
equalities. These also reduce to the binning results 
asymptotically as s = in the triple binning case and 
as g — >■ cxD or m -^ in the tanh case, as shown in Fig. |3l 
Small q is less optimal as shown by dashed curve. How- 
ever, for the MABK inequalities the triple binning case 
(s 7^ 0) is less optimal than the standard binned case 
(s = 0) even if we select a nonzero optimal value of s. 
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Optimal Function 




Figure 4: (Color online) Maximum violations of functional 
CV inequality with GHZ states as a function of the number 
of modes A'^. Violation of the Bell inequality is achieved when 
B > 1. The violations using the optimal function (solid) are 
much stronger than the CFRD result (dashed). 



C. Optimized functional CFRD 

In order to get the strongest violation of LHV theories, 
we can optimize the function of observables by using vari- 
ational calculus: 



6B 

Wn 



6B_ 

Sgn 



s 



tN 



= 0. 



(nr=i{/«w)'+5«w)'})j 

(5.9) 
For simplicity, we assume the functions /„ and (?„ are 
odd. The numerator can be maximized by choosing or- 
thogonal angles, while the denominator of the fraction is 
invariant with angles. In the case of the GHZ states, we 
find that: 



B 



N 



2^-^(f)-(nti/,t+n:=i4") 



n 



?i=i 



• 1 lfe=r+l -'n T^ 1 ln=l ^n 1 lk=r+l 



where I^ — 2 J e ^^ xf^dx, /„ 



= Ajx'^e^'^'"' 



In 

(5.10) 

m? + 



if-ndx, and /° = /e-2- [{f+f + (/-)2)]dx (given 
in the Appendix) are different integrals for x which con- 
tribute to the expectation values in both sides of inequal- 
ity (|5.8p . Here /^ = /„ ± 5„, and the factor e~^^ was 
obtained from the joint probability of observables. Re- 
quiring SB^/Sfn — 0, we find the optimal condition: 
fn{x) — ±(7„(x). The components of complex functions 
fm g-n are the same at each site, and have the form 



fn{x) = gn{x) 



1 + SnX^ 



(5.11) 



For the even TV case, it is optimal to choose r = N/2. 
Then em is independent of TV, but has to be calculated 



numerically since it satisfies a nonlinear integral equa- 
tion: En — 4:1"/!. Using the method shown in the 
Appendix, we can obtain the fixed values of integrals. 
We find that en = 2.9648 gives the optimal case, and 



B 



N 



■)N-2 



2(/+)V(^/°/) 



N/2 



For iV an odd number, the greatest violations occur 
for r = {N — l)/2. The optimal function has the same 
form as in (j5.1ip except that the parameter Epf changes 



to 



•AT, 



where: 



e'n = £n[ 



Ne+ 



-N^ 



Ne^ 



(5.12) 



-Af 



jY — En i 4. However, the numerical value of en 



and £^ 

and e^ now depend on N, as the integral equation (|5.12p 
for odd values of iV is A^-dependent. This provides bet- 
ter violations of (|5.8p than any other arbitrary function, 
provided TV > 5. The maximum B^ value with this opti- 
mal choice is shown in Fig. |4l compared with the CFRD 
result which uses a simple correlation function. 

In summary, a functional optimization approach gives 
a substantially larger violation of the CFRD inequality, 
and we have carried this optimization out explicitly for 
the case of a generalized GHZ state. 



VI. MANY OBSERVABLES PER SITE: SV 
INEQUALITIES 

In this Section, we compare the CFRD inequality with 
the inequalities given by Shchukin and Vogel (SV) |19J . 
These generalize the CFRD inequality to measurements 
with 4 or 8 settings or observables at each site. In princi- 
ple, this may allow larger violations or greater robustness 
against loss and noise. However, these are increasingly 
more complicated both to analyse and to carry out ex- 
periments for. By comparison, the CFRD approach con- 
siders just two detector settings per site. 



A. Four setting SV inequality 

First, we evaluate the results of SV for 4 observables 
per site, using arbitrary choices of phase for the quadra- 
ture measurements, where F„ — (A„")™ + i(A'„")™ -|- 

g3 gi 

j(X„")™ + fc(Ar„")™, are quaternionic functions of ob- 

gl 

servables. Here Xn"' , (i — 1, 2, 3, 4) are real ob- 
servables at site n. The multiplication rules for the 
quaternionic units i, j, and k are i^ = j^ = fc^ = —1, 
ij = —ji = k, jk = —kj = i, ki — ~ik = j. Thus, 
Shchukin and Vogel (SV) fl9'| obtain the following in- 
equality within any hidden-variable theory: 



\{FiF, 



2---/PI 



< 



(E(^i^^ 



1=1 



{x\ 



'2\2m 
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We consider the class of entangled states of type p.4p 
to find violations of this inequality. We use the optimal 
parameters found for the CFRD inequality, r = N /2 and 
ci = C2 = l/v2, for comparison with the SV inequalities 
with 4 observables per site. 



1. Two-site case 

As these inequalities are more complicated theoreti- 
cally, we begin with the simplest example of TV = 2. At 

site n = 1, we have 4 choices of observables [Xj^)"^ {I — 

1, 2, 3, 4) with corresponding outcomes I a; ^^ I ""sign (x^^), 
and similarly for the other site. The values of LHS and 
RHS are evaluated in terms of integrals for a;. First 
we evaluate the inequalities for A^ = 2, with the state 
IV') = ci|0)|l) + C2|l)|0), using ci = C2 = 1/V2 and 
m— 1/3. The LHS of the SV inequality has many terms: 



((X^')i(X^')i) = 0.42583 are obtained following the 
methods given in the Appendix. We find the value of 
the LHS of Shchukin and Vogel's inequality is 

LHS = I 0.434881{[cos((/)i-6ii)-cos((/)2-6'2) 
-cos{(l)^ - 6*3) - cos{(j)^ - e% 



COS[(p ~ 

-e^) + cos{(p^ 
- e-^) - cos{( 

q2\ I „„„/i,4 



0') 



4 a3\ 



+i[cos[i 
—cos{<j)'' 
+3[cos{ 
+cos{(j)'' 
+k[cos{ 
+cos{(j)" 

using angles like 9^ = 
-(j)3 = -(/,4 = _7r/2, to get LHS =\ 0.434881 x 4 p as aU 
imaginary values are 0. There are many optimal choice 
of angles, but the maximum values of the LHS are same. 
The value of the RHS is 



0, and e^ = 9^ = i 



(6.4) 



LHS 



\{[{X'^)' 


- + z(xf 


)i+jix2fs+kiXi') 


x[(^f) 


^ + i{xf 


f^+j{xf)^ + k{xf 


1 {[{xi) 


Hxi)^ 


-{xf 


Hxf)i 


-(^f)^ 


ixf)i - 


(^f) 


Hnr^] 


+z[{X<^2) 


Hxt)'^ 


+ {x<2 


)Hxtt^ 


-ix^)^ 


(^f)^4 


<x2) 


Hxt)'n 


+Mx2 


rHxtf^ 


+(^r 


)Hxt)'^ 


Hx'l)'^ 


ixf)^ - 


- ix2) 


Hxi)'n 


+M(^r 


)Hxt')^ 


-in 


)Hx^)i 



)m' 



+{x':)Hxff^+ix':)Hxt:f^])\' 



(6.2) 



RHS = 16((X^')2/3(xf )2/3) 
= 0.42583 X 16 , 



(6.5) 



so that the Bell ratio is the same as with the functional 
CFRD approach, using a power law with m = 1/3: 



B. 



2-SV(4) 



I 0.434881 X 4 |2 
0.42583 X 16 

I 0.434881 X 2 |2 
0.42583 X 4 

B2-CFRD ■ 



(6.6) 



This is caused by the optimal value of the phases cor- 
responding effectively to only two observables per site, 
since three of those phases have the same value. 



Similarly, the RHS is: 
RHS = ({(X^')2/3 + (xf)2/3 + (xf)2AV(xf)2A^} 

{(Xf )2/3 + (If )2/3 + (If )2/3 + (If )2/3}) 
= ((X«')2/3(xf )2/3 + (X^')2/3(xf )2/3 



^B 



Hx':f"{xtf/'+{x^:f/\xi; 

+ (Xf )2/3(xf )2/3 -f (If )2/3(xf 

Hx':f'Hxtf" + {xtf/\xt 
Hxtf/'{xt?" + {xtf"{xi 

+(xr)2/3(xf)V3+(xr)2/3(xf 

+{xtf/^{xtY" + {x<!:Y"{xt 



2/3 
2/3 
2/3 
2/3 
2/3 
2/3 

'") , (6.3) 



where we denote sites n ~ 1,2 as A^B respectively, 
with phase angles 9\ = 9^ and 62 — 4>^ ■ The av- 



erages of {{X'^^)i{xi)^ 



2. Four site case 

Next, we check mode number A^ = 4, and r — 2, for 
states like \Tp) = l/^/2(|0011) -t- |1100)). In QM, the 
operators at each site in Shchukin and Vogel's theory are 

Fa = (x«')"-f*(if)™-Hj(^rr + M^rr, 
fb = {xt^r+i{xfr+j{x'ir + k{xtr , 
Fc = {x^cT + ^{xi^r + jix^'r + Kxtr , 



fd = {XI r + Kx't, )" + 3{x% )'" + fc(xs r 



(6.7) 
7' , and 9'-^ = (pK After 



0.433881cos( 



9') and 



where again 9[ = 9\ 9^^ = (t^\ &3 
extensive algebra, we find that the largest Bell values 
with optimal choice of angles are still the same as in 
the CFRD case using m = 1/3: B4_sv{4) = B^-cfrd- 
The reason for this identity, as in the two-site case, is 
that with optimal settings we have effectively only two 
measurements per site. 
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B. Eight setting example 

Finally, we can check the SV inequality with 8 ob- 
servables at each site. This approach uses octonions, 
and results in a rather complex inequality |l9| . There- 
fore, we only consider here the simplest two-site case 

IV) = Ci|0)|l) + C2|l)|0), Ci = C2 = I/V2. 

We consider the fractional case m = 1/3 as an example. 
The LHS of the inequality is: 



|O.434881{[cos(0i-6|i)- 


- cos{(p'^ 


~o') 


-cos{(l)^ - 6*3) ~ cos{<j)^ - 


9^) 




~cos{(tf' - 9^) - cos{(j)^ - 


9^) 




-cos{(f ~ 9^) - cos{(j)^ - 


0')] 




+ii[cos{(f)^ - 9"^) + cos{4>'' 


'■-9') 




+cos{(j}^ -9^)- cos{(j)^ - 


9^) 




+cos{4>^ - 9^) - cos{(l)^ - 


9^) 




-cos{4>' - 9^) -f cos{cf)^ - 


9')] 




+i2[---]+h[---]+ii[---' 


\+i5[--- 


] 


+^,[■■■]+^l[■ ■■]}?, 




(6.8) 



where Fa = (X^ )"'-t-ii(X| )'^+i2{X\^ )'^+h{X'{ )"-|- 
i4(^f)"-t-i5(xf)"+i6(^A)™ + »7(^f)'",isanocto- 
nionic function of real observables, (X^ )™, (Z = 1 — 8) 

at site A; similarly for site B, involving (X^ )'". The 
multiplication rules for the imaginary units of octonions 
were listed in [19]. One of the optimal choices of an- 
gles still keeps the same pattern as in two-sites case: 
6»i = 01 = 0, and 9^^^ == -^'^^ = -7r/2, so that we 
obtain LHS =\ 0.434881 x 8 p, and 



RHS = 64((X^')2/3(lf )2/3) 



Therefore 



B 



2-51^(8) 



|0. 434881x81^ 
0.4258x64 



B2-CFRD 



(6.9) 



(6.10) 



strategy that may allow genuine loophole-free violations 
of LR. However, it may be argued that since increasing 
the number of sites will increase the number of detectors 
required, there will be no advantage. Only careful calcu- 
lation of the Bell observable Bm including the detection 
efficiency ry can determine whether the strategy is advan- 
tageous for loophole-free violation of Bell inequalities. 



A. Model of loss and impurity 

Loss is modelled using a beam splitter as follows. The 
field modes ai„ at each site are independently coupled 
to a second mode a^,, assumed to be in a vacuum state. 
Bosons are lost from the field into the vacuum mode, the 
strength of coupling determining the rate of loss. For 
each mode, the beam splitter model gives the final de- 
tected and vacuum mode in terms of the inputs ai„ and 
a,. 






(7.1) 



where 77 is the efficiency, the probability for detecting a 
boson, given that one is incident. One can rewrite these 
relations to obtain 



= y/rjaout + \/l - r]0'v,out 



(7.2) 



Since we only measure "dout", not ^^Uv^out", we need to 
trace over all the unmeasured vacuum output modes to 
obtain the final density operator for the detected modes 
after loss. 

Thus, given a total density matrix of /5, we only mea- 
sure a reduced density matrix pout, where: 



Pout — Tr„^out{p} 

= y^^('ipk\v,outp\lpk)v,o 



(7.3) 



For A^ = 4 modes, the optimal result of the SV inequal- 
ity with 8 observables at each site is also the same as for 
the functional CFRD fractional moment inequality, for 
the particular state \^|J) = 1/V2(|0011) + |1100)). 

In summary, in the cases investigated, the results are 
the same as those found with the simpler CFRD ap- 
proach, using fractional moments. The reason is that the 
optimal values can be found with only 2 effective settings 
per site with the types of states we consider. However, 
this situation may change with more general input states 
or functional transformations. We leave this as an open 
question for future work. 



We also examine the effect of impurity, by considering 
an input state p' = p\'ijj){'ip\ + (1 -p)pmix, where p^ix is 
the mixed state obtained by complete decoherence in the 
occupation-number basis, i.e. 

Prmx = |cinO)«ni>^^'^"'^^(0|^''(l|^^'^'''^ + 

+ |c2p|l)®''|0)®(^-'')(l|^''(0|®'(^-'') . (7.4) 

Here p is the probability the system is in the original pure 
state p.4p . While this is not the most general model of 
state impurity, it is a relatively simple one that allows 
us to carry out indicative calculations of the effects of 
impure input states. 



VII. SENSITIVITY TO LOSS AND MIXING 



B. Three- mode MABK example 



The value of the Bell observable B^q or S^ increases 
with the number of sites N , and this is suggestive of a 



We can transform the input state for A'' = 3 and r = 1 , 
IV) = ^(|100) -f- |011))|000), for example, where |000) 
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represents the vacuum fields, into the output after loss. 
We denote the three field mode operators involved as a, 
6, c for clarity. Thus in this case, 

\out) = [aL|000)|000) + 5^4 1000) |000)]/V2 
= (V^|100)o„t|000),,o„t + 

^/l-V\000)out\lOO)v^out + ??|011)o«t|000).,o„t 

+(i-v)\ooo)out\on)y^out + 



^/T^^i\00l)out\0l0)v,OUt 

+ \01Q)out\001)v^out)/V2 . 



(7.5) 



Next, we need to trace over the vacuum output modes to 
obtain a reduced density matrix pout, where: 

Pout = Tr„_o„t{/5} 

= ^^{ijk\y^outp\ijk)v,out 

= ^{[V^llOO) + r;|Oll)][0y(lOO| + 77(011|] + 

(1 - ?7)|000)(000| + 77(1 - r;)|001)(001| 

+77(l-7y)|010)(010| + 

+(1-77)2|000)(000|} . (7.6) 

Proceeding, the probability distributions are evaluated. 
Thus, in this three-mode example: 



P{xe,x^,Xj) = {xe\{x^\{xy\pout\xe}\x4,)\x^} 

= i(l)3/2g-2(3;2+2:J+2:2) ^ 
2 TT 

X [ I {2^e'''xe+ive"^e'^x^x^) p + 

+47j{l - Tj) X { {\ x^ \^ + \ x^ \^} 

+ [l-7/+(l-r/)2]] , (7.7) 

and the integrals give for example 



>0,bin -^(p,bin -y'7,6'in\ ( ■^'7 



{XY'^X^'^'^Xy'") = (^j cos(7 + 0-0)(7.. 

Examining the MABK Bell inequalities defined by 
(|4.2p . we finally obtain, using the "orthogonal" choice of 
measurements given for (|4.7p . the result: 



277 

I {oliJjbinCbin) I = 4 I — 



(7.9) 



Alternatively, for the measurement angles 9 — Q, (j) — 
_7r/6, 7 = -27r/6, 9' = 7r/2, 0' = -47r/6, 7' = -57r/6, 
we can rotate the LHS moment, but so that again 



'At 



alinbbinCbin) H 4 
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(7.10) 



C. Higher N results for MABK inequalities 

Continuing in this fashion for higher N, we find that 
the values of the Bell observable including the effect of 
detection inefficiencies with the optimal choice of angles 

is 



S, 



N 



1 U-q^ 
V2 V~ 
which implies a critical minimum efficiency 



(7.11) 



(7.12) 



in order to violate the MABK inequalities (|4.2p . 

We can also examine the effect of impurity or noise, 
by considering a state p' — p|7/;)('0| + (1 — p)pmix, as 
described above. The effect of this type of 'noise' on the 
MABK inequality is to give 



SNip,v) 



V2 fAT]p^ 



(7.13) 



This result agrees with that obtained by Acin et al [T^ 
for the case r] — 1. 

For lower N, the strategy of binning and using the 
MABK inequalities shows an advantage, by allowing a 
violation of LHV theories for iV = 3, 4, 5. However, 
very high efficiencies, r] > 0.99, 0.93, 0.90 respectively, 
even for a pure state p = 1, are required, as shown in 
Fig. [5l While high detection efficiencies are feasible for 
homodyne detection, these efficiency values are still quite 
challenging once generation losses are also taken into ac- 
count. In view of this, the high requirement for r]mm for 
the case A^ = 3 may be prohibitive. 



0.08 



0.06 




0.04- 



0.02 



0.9 ri . 0.95 
'mm 



Figure 5: (Color online) The boundary of maximum noise 
(1 — p)max and minimum detection efficiency rjmin to violate 
the MABK Bell inequality with binned outcomes. Different 
TV is distinguished by the line style. 
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efficiency. The violation is less sensitive to detector inef- 
ficiency in the macroscopic, large N limit. 

Compared with the MABK binning approach, the ro- 
bustness of the functional inequality against noise and 
inefficiency is apparent as shown in Fig. [S) For iV > 7, 
the functional CV inequalities used with an optimal func- 
tion shows clear advantages, allowing violation of LHV at 
much lower efficiencies and larger maximum noise. For 
example, for N — 40, we can get Bell violation with 
'^cr^t ^ 70.8% for pure state, or Pmax ^ 0.25 for ideal de- 
tection efficiency. The potentially feasible r\crit ^ 80% 
requires iV sa 10 for a pure state if the functional op- 
timization is used. By comparison, the binned case re- 
quires iV ^ 40, and the original CFRD cannot reach this 
even with N ^ oo. 



Figure 6: (Color online) The boundary of maximum noise 
(1 — p)max and minimum detection efficiency rimin to violate 
the functional CFRD Bell inequality with observables as the 
optimal function. Different A*' is distinguished by the line 
style. 



D. Loss and noise with functional inequalities 

This approach is also applied to enable a prediction of 
the effect of loss and noise on the functional inequalities 
of Section IV, and the results are plotted in Fig. IHlfor the 
optimum function. Here, the parameter sn is changed to 
SNiv) ior the optimum function. For N even we find that 



Bn 



2-qeN 



27] + (l-r;)ew ' 

Ml+)\r^p)\. 



^N~2\ 



7r/°C 



(7.14) 



where en is defined as before, and C = riI-\-{\ — rj)I° . For 
the case of odd N the relevant integral equations change, 
giving a modified (and slightly reduced) Bell variable i?^, 
where: 



2\/7°C 



B' 



N 



I" + C' 



:B 



N 



(7.15) 



Here e^{T]) = SNiv) + 4, and Bn is defined as in Eq 

For pure state, in the case of r = N/2, N is even, we 
can give the threshold efficiency requirement analytically 
since all integrals are independent of TV, where 



_ [(4-£Ar) + V(4-£^)2+4£^]/o/^ 
Vent- 26-4/W(-/+-)4 ■ V-^"-* 

This reduces at large N to an asymptotic value of rj^ = 
0.69, thus dramatically reducing the required detection 



VIII. CONCLUSION 

In summary, we have developed a functional approach 
to tests of quantum nonlocality, and showed how this 
gives an exceptional degree of robustness against losses 
and noise. For continuous variable measurements on 
GHZ states, this approach gives distinct advantages when 
compared with either binning methods or the moment- 
based approach of Cavalcanti et al. (CFRD) [18|. We 
have also compared the class of inequalities derived by 
CFRD with the generalizations given by Shchukin and 
Vogel (SV) for more observables at each site but found 
no better violations with the photon-number correlated 
states introduced here. 

However, there are still many open questions. The first 
is applicable even to the CFRD result: we do not know 
yet if there are any other quantum states that can give 
better violations than those given here for the GHZ-type 
states. Another question is whether there are other states 
for which the multiple-observable SV inequalities provide 
an advantage over the CFRD. One may also ask if there 
are yet more general continuous variable inequalities than 
either the CFRD or SV resufis. 

Further applications of these ideas may include appli- 
cation of the functional approach to tests of other forms 
of nonlocality — e.g, entanglement and EPR-stccring 
|3a |. It seems possible that a functional optimization 
will allow more robust implementations of these tests as 
well. 
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Appendix I: Quadrature Wavefunctions 



B. General rotation 



Calculations with the functional transformations used 
here rely on the fact that quadrature variables have the 
same wave-functions as the quantum harmonic oscillator. 
We first recall that the harmonic oscillator Hamiltonian 
with unit mass and frequency lo can be written: 



H^\{P^+uj^x^) 



.1) 



Here p is the momentum and x is the position operator, 
with commutators [x,p\ = ih, and creation and annihila- 
tion operators defined as: 






(8.2) 



For simplicity we assume that u) = 1 and h = 1/2, which 
corresponds to the convention used throughout the pa- 
per: 

X ^ {a) + a) 12 , 

p = -i{a-a^)/2 . (8.3) 

A. Wavefunction for quadrature variables 

In this standard form, the harmonic oscillator wave- 
function |3J| for the m}^ energy eigenstate is: 



{x\n) = ( 



7r22"-i [n!] 



;f'^H^{V2x)e-^" 



(8.4) 



where Hn{x) is the n*'' Hermite polynomial. Thus, for 
example: 



(x|0) 



(x\l) 



1/4 



1/4 

2a; ( - ) e 

TT. 



(8.5) 



We also need the complementary (momentum) wave- 
functions: 



(pIV) 



1 



exp{~2ip'x'){x'\ii}dx' , (8.6) 



Again, using normalized quadrature variables, one finds 
that: 



(P\0) 



1/4 



{p\l) = -2ip ( - I e 



1/4 



(8.7) 



Next, we consider the case of a general rotation of the 
quadrature phase by defining 

X^ = xcose+psme = {ae-'" + a^e:'^)/2 , 

Hence we find that 



(^ |0) = Z 



(A'^ll) = 2e 



2\ ^_^xOf 



1/4 



xSp-ix')" 



and p corresponds to 9 = tx 12. 



C. Joint probabilities 

To carry out moment calculations, it is necessary to 
have joint probabilities. As an example of this calcula- 
tion, we consider a two- mode state p = |^) (^|, where: 



|vl/) = i=[|0,l)-t-|l,0)] 



(8.9) 



Carrying out the calculation of the joint probability, 
and defining |X|^ = {X^f + {X'^f, we find that: 

P(X\X'^) = {X'^\{X'f'\p\X^)\X'^) 



-|(A«|0)(A^|l) + (A«|l)(A^|0)p 

-)|e-l^l'(e^^A* + e^''A'')|2 
n J 



1 ^-2|X|^ 



]X\^ + 2X'^X'''cos{0-cl))] (8.10) 



D. Binned moment calculation 

As an example of the application of the joint proba- 
bility, consider a binned moment calculation for the case 
N = 2, using an initial Bell state of the form given above. 
Then we obtain: 






OO /'OO 



p{x\x*)d^y;. 



X0=0 JX'I>=Q 
OO M 



/•OO />u 

/ / P{X\X'^)<fX 

Jxo=OJX'l>=-oo 

/ / P{X\X'f')<fX 

/ / P(X^,X''')(fX 

(8.11) 



- cos(0 - 9) . 

TT ' 
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mized by using orthogonal angles so that cos{ 



±1: 



E. Functional moment calculation 

Here we show how to calculate the expectation value 
of products of functional moments using the joint prob- 
abilities. For the case N = 2, using the same state given 
above, the value of the LHS is 



LHS 



f[[UX:) + ^g,,{X'J] 



\n=l 



.12) 



This can be evaluated by the joint probability shown in 
Eq. (|8.10p . and will include terms like: 

{hixt)Mm = (1) J J e-'i^n,M\xf 

+ 2X''X'f'cos{d - 0)]<i2x . (8.13) 

To show this more clearly, we give an example of the 
observables as fractional order moments (X^)™ with cor- 
responding outcomes \x^\"^sign{x^). Replacing /i and /2 
in the above integral with (X^)'" and (X^)™, we obtain 
that 



((11)^(1^)3) 

{ix%)Hxt)h 



0.433881cos( 
0.42583 



.14) 



for m = 1/3 and m ~ 2/3. 

For convenience, we assume the functions /„ and gn 
are odd, and thereby only the cross term in (|8.13p has a 
nonzero contribution. The value of LHS can be maxi- 



LHS 



nL l(/" +9n)+ Ul=lifn - 9n) 

2 2 

n=l 



n— 1 n— 1 



15) 



where f^ — fn ± g„, and define integrals for x: In = 
2je-^-\f^dx. 
While the value of RHS 



:n [/"(^n)' 



9'\2 



RHS = (II [fnixy+gUX^^r ) 

n=l 



/TT ifn + 9n) + [fn ~ 9n) ^ ,„ ^ „s 



ri=l 



is invariant with angles as it only related to the square 
items in the probability of Eq. (|8.10p . Define /„ = 
4/x2e-2-^[(/+)2 + (/-)2)]dx,and/° = /e-2-^[(/+)2 + 
{fn )'^)]d'X, then we can obtain 



^ r N r N 



n=l k=r+l ri=l k=r+l 



(8.17) 



where r = 1, iV = 2. 
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